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Abstract 

We demonstrate that a direct approach to improving Hamiltonian lattice gauge 
theory is possible. Our approach is to correct errors in the Kogut-Susskind 
Hamiltonian by incorporating additional gauge invariant terms. The coeffi- 
cients of these terms are chosen so that the order a 2 classical errors vanish. 
We conclude with a brief discussion of tadpole improvement in Hamiltonian 
lattice gauge theory. 



1 Introduction 

The idea of using the lattice as an ultra-violet regulator for quantum chromodynam- 
ics (QCD) was proposed by Wilson in 1974 in his action formulation of lattice gauge 
theory [jij. Soon after, Kogut and Susskind formulated the corresponding Hamiltonian 
version of lattice gauge theory ||]. Both approaches were developed by demanding the 
correct continuum limit be obtained in the limit of vanishing lattice spacing a. Creutz 
showed that the Kogut-Susskind Hamiltonian could be derived from the Wilson action 
using the transfer matrix method ||. Later, Kogut demonstrated that the same could be 
done by taking the continuous time limit of the Wilson action and performing a canonical 
Legendre transformation 

To date, the majority of work in lattice QCD has been performed in the action formulation. 
An advantage of this approach is that it readily lends itself to Monte Carlo techniques. 
Working in the Hamiltonian approach brings a different intuition to the problem and 
serves as a check of universality. An advantage of Hamiltonian lattice gauge theory is 
in the applicability of techniques from many body physics Q. Also, it appears that in 
finite density QCD, a Hamiltonian approach is favourable due to the so-called complex 
action problem which rules out the use of standard Monte Carlo techniques in the action 
formulation Q. 
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Much work in the past decade has been devoted to improving lattice actions 0. Initiated 
by Symanzik in 1983 @, the aim of the improvement programme is to reduce the deviation 
between lattice and continuum QCD. For pure gauge theory on a lattice, the deviations 
between continuum and lattice theories start at order a 2 . The motivation for improvement 
lies in the fact that the computational cost of a lattice QCD simulation is proportional to 
a~ k , where 6 < k < 7. It is by far more efficient to build an improved theory than it is to 
work on finer lattices. The improvement programme has allowed accurate calculations to 
be performed on relatively coarse lattices and brought the most complicated calculations 
within the reach of today's most powerful computers. 

In contrast, the improvement of lattice Hamiltonians has only recently begun. Perhaps 
the most extensive treatment to date is due to Luo, Guo, Kroger and Schiitter Q who 
discussed the improvement of Hamiltonian lattice gauge theory for gluons. In their study 
it was discovered that deriving an improved Hamiltonian from a Symanzik improved ac- 
tion, whether by transfer matrix or canonical Legendre transformation, results in a kinetic 
Hamiltonian with an infinite number of terms coupling lattice sites which are arbitrarily 
far apart. To derive a local kinetic Hamiltonian coupling only nearest neighbour lattice 
sites it was found necessary to start with an improved action with an infinite number of 
terms, coupling distant lattice sites. 

With this technique the order a 2 errors are removed from the Kogut-Susskind Hamiltonian. 
However, generating Hamiltonians with further improvement would seem exceedingly dif- 
ficult. This is because one would need to start from a Liischer-Weisz improved action 



with non-planar terms [10]. For this reason we propose a move to the Symanzik approach, 
as applied to the Hamiltonian. That is, to construct improved Hamiltonians directly by 
adding appropriate gauge invariant terms and fixing their coefficients so that errors are 
cancelled. 

To date we have implemented Symanzik improvement to order a 2 in the pure lattice gauge 
theory Hamiltonian (and to order a 4 in the kinetic part). We report those results here^]. 



2 Errors in Lattice Gauge Theory 

Before discussing Hamiltonian improvement we must first understand how deviations be- 
tween lattice gauge theory and its continuum counterpart arise. The deviations can be 
separated into two classes, classical and quantum errors, which will be described in what 
follows. 

Rather than being constructed from gluon fields, pure gauge theory on the lattice is built 
from link operators 

U^{x) = exp lig / dx ■ A \ . (1) 

I ./Link J 

3 The 0{a 2 ) results were given in preliminary form in 
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Figure 1: The diagrammatic representation of the link operator U^{x). 

Here g is the QCD coupling, [i = 0, . . . , 3 is the Dirac index labelling the time-like and 
space-like directions, x labels a lattice site and A is the gluon field. The integral runs 
along the link joining the lattice sites x and x + afj,. This leads to the diagrammatic 
representation of the link operator shown in figure |l[ Lattice gauge theories are built from 
link operators rather than gluon fields in order to maintain manifest gauge invariance. 

We define the lattice gluon field Aj^(x') to be the average of the continuum gluon field A 
along the link joining x and x + a^i: 

Afa') = - ! dx-A U„(x) = e i9aA ^ x '\ (2) 

a JLink 

where x' is a point near the points x and x + a^i. On the lattice the gluon field is defined at 
only one point along (or nearby) a link. This leads to interpolation errors in the integral 
in eqn For instance, by choosing to evaluate the gluon field at the midpoint of the link, 
the lattice and continuum fields are related by 

= a m (x) + yS Mx) + 4rS Mx) + • • • (3) 

We see that the lattice gluon field reduces to its continuum counterpart in the continuum 
limit (a — ► 0), but that they differ by terms of order a 2 , which may be called interpolation 
errors. 

Having discussed classical errors we now move on to quantum errors in lattice gauge theory. 
Quantum errors arise in two different contexts. Firstly, the lattice acts as an ultraviolet 
regulator allowing the simulation of only those states with momenta less than it/a. The 
absence of high momentum states results in a deviation between lattice and continuum 
theories. Secondly, non-physical interactions arise due to the use of the link operator in 
constructing the lattice theory. To demonstrate this we expand the link operator in powers 
of g, 

a 2 a 2 

U„(x) = 1 + igaA^x) - ^-A^A^x) + ..., (4) 

and note that the interaction of any number of gluons is allowed. Naively, the unphysical 
interactions are suppressed by powers of a. However, when contracted, products of pairs 
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of gluon fields produce ultraviolet divergences (oc f/a 2 ) which exactly cancel the a depen- 
dence of the expansion. These terms can be uncomfortably large and result in what are 
known as tadpole errors. 

In the last decade the improvement programme has led to a good understanding of both 
classical and quantum errors in quark and gluon actions (See reference || and refer- 
ences within). In contrast, only the lowest order classical errors have been corrected in 
the Kogut-Susskind Hamiltonian. Conjectures have been made about the structure of a 
quantum improved Hamiltonian ||, but a perturbative study has not yet been carried out. 

We briefly describe how quantum improved Hamiltonian can be constructed in § ||. 



3 Symanzik Improvement of the Lattice Hamiltonian 

3.1 Introduction 

In this section we derive an improved Hamiltonian directly using the Symanzik approach 
of adding irrelevant terms and fixing their coefficients in order to cancel errors. As a first 
step we aim to correct the classical order a 2 errors arising in the lattice Hamiltonian for 
pure SU(iV) gauge theory. 

The Kogut-Susskind Hamiltonian for pure SU(iV) gauge theory on the lattice is given by 

H®=kW+V®, (5) 
where the kinetic and potential terms are given respectively by, 

K (0) = ^Tr{£f(*)i^(*)} (6) 



2 
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Here E L is the lattice chromo-electric field, N is the dimension of the gauge group, and 
Pij(x) is the plaquette operator in the plane, 

^•(*) = l-^ReTr jQjj. ( 8 ) 

i 

We discuss the improvement of the kinetic and potential terms separately in what follows. 



3.2 Improving the Potential Term 

To improve the potential part of the Kogut-Susskind Hamiltonian, we follow the process 
of improving the Wilson action. Our first step is to expand the potential term in powers 
of a. This is done using Stokes' theorem to expand the plaquette operator, 

Pij(x) = 1 — — ReTr exp i.ig <j> A ■ dx 
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Tr {F^F^x)} + ^Tr {F tJ {x){D 2 + D 2 )F l3 (x)} + 



0) 



Here Di = di — igAi is the gauge covariant derivative, and are the spatial component 
of the gluon field tensor. The potential term then becomes: 



2N 
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x,i<j 



= a 3 E h^j(x)F i3 (x) + yTri^(x)(A 2 + D 2 )F i3 (x) + ... 

x,i<j ^ 

« ^ |^X>{i^(*^ (10) 

We see that the correct continuum limit is restored in the limit as a — > 0. To cancel the 
order a 2 error we introduce rectangular loops into the potential term. These are the next 
most complicated gauge invariant constructions on the lattice after plaquettes. We write 
the improved potential term as the linear combination 



ag 2 ^ 

x,K] 



Y 

XP i:j (x) + - (Rij(x) + Rji{x)) 



(11) 



where Rij (x) is the rectangle operator in the plane with the long side in the i direction, 
given by: 



%(x) = 1 - ^ReTr 



(12) 



The constants X and Y are to be fixed so that the order a 2 errors vanish. 



The rectangle operator can be expanded in powers of a using Stokes' theorem: 
Rij(x) ■■ 



4« 2 n 4 4-n 2 n 6 

-|^Tr {F i3 {x)F %J {x)} + ^Tr {ffc(s)(4D? + D^F^x)} + ... (13) 



Substituting eqns |9| and into eqn 11 we see that cancelling the order a 2 error in 
requires X = 5/3 and Y = —1/6. This leads to the improved potential term 
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(i^-(x) +-R,t(a?)) 



(14) 



In principle, the next lowest order classical errors could be corrected by including addi- 
tional, more complicated terms in the potential term. This has not been done because 
many additional diagrams are required to cancel the large number of order a 4 error terms. 
Since these errors are swamped by order a 2 g 2 quantum errors, addressing quantum cor- 
rections in the Hamiltonian approach would seem to be of more immediate importance. 
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3.3 Improving the Kinetic Term 



Constructing a kinetic Hamiltonian with a finite number of terms has proven to be a 
nontrivial exercise. Luo, Guo, Kroger and Schiitter demonstrated an interesting trade off 
when using either the transfer matrix or Legendre transformation methods to derive an 
improved Hamiltonian Q. Both techniques require the starting point to be an improved 
action. When one starts from an improved action incorporating rectangular terms the 
resulting Hamiltonian has infinitely many terms and couples links which are arbitrarily 
far apart. To produce a Hamiltonian which couples only nearest neighbour links, it was 
found necessary to start from a carefully constructed highly non-local improved action. 

Here we demonstrate an alternative approach, similar in nature to the Symanzik improve- 
ment of the Wilson action. One only needs to include additional gauge invariant terms 
with appropriate continuum behaviour in the kinetic Hamiltonian. The coefficients of the 
additional terms are chosen so that the order a 2 errors vanish. 



An important step in understanding the errors that arise in the kinetic Kogut-Susskind 
Hamiltonian involves making the distinction between lattice and continuum fields. In § |2| 
we defined the lattice gluon field on a link to be the average of the continuum gluon field 
along the link: 

Afe) = \J dx-A = A,(x) + ^flA^x) + ^d*A^x) + . . . (15) 

From this we can build the following sequence of approximations to the lattice gluon field: 

A?\x) = A t (x) 

Af\x) = Ai (x) + ±a 2 d?Mx) (16) 
Af\x) = Al (x) + ^a 2 d 2 Ai{x) + -^a'dtMx) 



Perhaps the most important property of the electric field is that it generate group trans- 
formations. Mathematically, this translates to the electric and gluon fields satisfying the 
commutation relations, 

[E?(x),AP(y)] = -^6 xy 6 ij 5 a p. (17) 

It is desirable for this hold on the lattice for any degree of approximation. Let us consider 
what happens to these commutation relations on the lattice for the approximation labelled 



by the superscript (1) in eqns 16: 



[E^ a {x),Af\y)] = [E^ a (x),A^(y) + ^flA%)\. (18) 

We observe that if the lattice electric field is taken to be the continuum electric field, order 
a 2 errors arise in the commutation relations. To cancel this error we set 

E^{x)=E?(x)- ( ^d 2 E?(x). (19) 
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We can take this to order a 4 by setting 

E^ a (x) = E?(x) - ^d 2 Ef(x) + J^dfEfix). (20) 
In this way a sequence of approximations to the lattice electric field E L can be constructed. 

Making use of these approximations we can analyse the classical errors arising in the ki- 
netic Hamiltonian. To cancel these errors we take the approach of adding new terms and 
fixing their coefficients in order to cancel the order a 2 error. We have a great deal of 
freedom in choosing additional terms. They are restricted only by gauge invariance and 
the need for an appropriate continuum limit. 

To understand the construction of gauge invariant kinetic terms it is important to recall 
that the electric field and link operator transform as follows under a local gauge transfor- 
mation A(x): 

Ei{x) -» A(x)Ei(x)A\x) (21) 
Ui(x) -> k{x)Ui{x)tf(x + ai). (22) 

Consequently, the next most complicated gauge invariant term we can construct (after 
Tr E L E L ) couples nearest neighbour electric fields: 

Tr {Ei(x)Ui(x)Ei(x + ai)u}(x)} . (23) 

More complicated gauge invariant terms are easily constructed. One only needs to couple 
electric fields on different links anywhere around a closed loop. Consequently, generat- 
ing Hamiltonians with higher degrees of improvement would be seem to be more readily 
achieved within this approach. 

Incorporating nearest neighbour interactions leads to the simplest improved kinetic Hamil- 
tonian: 

K m = ^Y,^{ XE i( x ) E H^+YEi'(x)U i (x)E('(x + ai)uj(x)}. (24) 

x,i 

We fix the coefficients X and Y to cancel the order a 2 error. To do this we expand the 
second term in a Taylor series in a . Ignoring 0(g 2 a 2 ) errors, we then substitute E L E^ 
from eqn To cancel the order a 2 error we must set X = 5/6 and Y = 1/6. This results 
in the order a 2 improved kinetic Hamiltonian 

Kil) = yE Tr {^(*) + \E?(z)U t {x)E?(x + ai)u}(x)\ . (25) 

x,i 

This is the result of Luo, Guo, Kroger and Schiitter Q. We can take the this to order a 4 
by including next nearest neighbour interactions. A similar calculation using E L « 
from eqn |2TJ gives 

x,i 

-^E[(x)Ui{x)Ui(x + ai)E^(x + 2ai)uj(x + ai)U}(x)\ . (26) 
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4 Tadpole Improvement 

Tadpole improvement, developed by Lepage and Mackenzie fi"2|| , is an important step in 
removing errors from lattice gauge theory. It is necessary for close agreement between 
lattice perturbation theory and Monte Carlo calculations on coarse lattices. 

In the action formulation tadpole improvement is handled by dividing all link operators 
by the mean link no- In the Hamiltonian approach two conflicting implementations have 
been suggested. The earliest starts from a tadpole improved action and carries factors of 
uq into the Hamiltonian |J. More recently it was suggested that no tadpole improvement 
was necessary in the kinetic term of the improved Hamiltonian [ |13| . Here we present our 
own views on the correct implementation. 

In the Hamiltonian approach the question of whether the electric field should be scaled 
arises. This question is easily answered by considering the commutation relations between 
the link operator and electric field: 

[E?{x), Uj(y)\ = ^SijS^Uiix). (27) 

We see that if we divide all link operators by uq we have 

[E?(x), —U 3 (y)} = ^Si^y^—Uiix). (28) 

Uq a Uq 

We observe that the electric field cannot be rescaled and still maintain the correct com- 
mutation relations. Thus under tadpole improvement the electric field cannot change. We 
must, however, divide the second of the kinetic terms by a factor of Uq. Tadpoles arise in 
this term because the electric and gluon fields do not commute. 



Including tadpole improvement in eqns 14 and 25 leads to the order a 2 tadpole improved 
Hamiltonian: 

" Tr^Et(x)Et(x) + -^EtidUiWEtix + ai)U}(a 



X,l 

2N 
ag 2 



E 

x,i<j 



— jPa(x) 7? (Rii(x) + RjAx)) 



(29) 



5 Conclusion 

We have demonstrated that direct improvement of the Kogut-Susskind Hamiltonian by 
demanding the correct continuum limit is possible. The advantage of our direct approach 
is that it is easily extended to more complicated Hamiltonians. One simply needs to 
construct suitable gauge invariant terms to add to the kinetic Hamiltonian and fix the 
coefficients so that higher order errors are cancelled. 



Our next step is to perform variational and coupled cluster SU(3) calculations to determine 
precisely the level of improvement achieved by the improved Hamiltonians. Other groups 
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have made progress with these calculations for U(l) [|l3| and SU(2) |l4j] with promising 
results. 



In the near future we intend to extend the direct approach to the cancellation of quantum 
errors which have not yet been examined in Hamiltonian lattice gauge theory. 
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